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In the present paper we initiate the study of the product higher rank numerical range. 
The latter, being a variant of the higher rank numerical range [M.-D. Choi et al, Rep. 
Math. Phys. 58, 77 (2006); Lin. Alg. Appl. 418, 828 (2006)], is a natural tool for 
studying construction of quantum error correction codes for multiple access channels. We 
review properties of this set and relate it to other numerical ranges, which were recently 
introduced in the literature. Further, the concept is applied to the construction of codes 
for bi— unitary two— access channels with a hermitian noise model. Analytical techniques 
for both outerbounding the product higher rank numerical range and determining its 
exact shape are developed for this case. Finally, the reverse problem of constructing a 
noise model for a given product range is considered. 
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1 Introduction 

Quantum information transmission [TJ [21 121 HI [5] inevitably involves occurrence of errors 
which faced not properly can disable faithful communication. Developing fruitful and useful 
strategies for combating these errors is thus one of the main challenges of the theory of 
quantum channels. Many effort has been put in this direction and several techniques have 
been developed to overcome destructive influence of coupling to the environment (see, e.g., 
E]). Among them, quantum error correction codes (QECC) along with the celebrated 
Kmll-Lanamme (KL) conditions [5] are most widely recognized (see [3] and references therein). 
Methods of constructing QECC for quantum communication have been previously reported 
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in the literature [TOj [ITJ [12]. There has also been a significant progress in experimental 
realizations of the propositions [HI HU [15] . 

Recently, a fruitful approach to solving Knill-Laflamme conditions has been put forward 
[16] [17] [18]. It relies on the concept of the higher rank numerical range of an operator 
and provides a systematic framework for a construction of quantum error correction codes 
(QECC). As shown in Ref. [El [17] it may capture more possibilities than the stabilizer 
formalism [12] , 

However, all proposals analyzed so far concerned bipartite communication (see however 
the related paper [21]). No general approach has been developed to treat the case of the 
larger number of users of a quantum network. The main purpose of this paper is to provide 
a mathematical tool of the product higher rank numerial range for designing QECC for mul- 
tiaccess quantum communication. Product higher rank numerical range by the definition is 
a higher rank numerical range restricted to product projections. 

The paper is organized as follows. First, we review basic ideas, including error correction, 
of quantum communication over (multiparty) quantum channels. Further, we briefly recall 
some basic features of the higher rank numerical range with applications to error correction. 
We then move to the main body of the present paper by introducing the notion of the product 
higher rank numerical range and discussing its basic properties. In the next section we 
consider construction of QECC for bi-unitary channels with a hermitian noise model. We 
then demonstrate exemplary applications of our findings to some concrete problems. We also 
discuss the reverse problem of constructing a noise model for which a product code exists. 
The manuscript is concluded with a discussion. 

2 Quantum channels and quantum error correction 

Here we briefly recall some basic ideas of communication over quantum channels and set the 
scenario for further considerations. 

2.1 Quantum channels 

Quantum channel £ is a completely positive trace-preserving map. Every channel admits 
the so-called Kraus (or operator-sum) representation as follows C(g) = Yli^iQ^H with 

A\Ai = 1 [22] [23]. A random unitary channel is the one which has the representation 
C(g) — ^iPiUiQU} , where Ui are unitary and YliPi = 1' Pi — 0- When such a channel has 
two Kraus operators, i.e., C(g) — pUigU^ + (1 — p)UigUi , it is called a bi-unitary channel 
(BUC). This kind of channels is the main interest of the present paper. 

Channels can be classified upon the number of senders and receivers using them. We 
have the following types of channels according to such a classification [4j [24] [25] [26l 127] : 
(i) bipartite — one sender and a single receiver, (ii) multiple access — several senders and 
one receiver, (iii) broadcast — one sender and several receivers, (iv) fcm-user — k senders 
transmit information to m receivers (fc,m > 1). 

In our reasonings we mainly concentrate on two-access channels, that is multiple access 
channels with two senders. 

Due to the possibility of a global rotation u\L{g)U\ on the output of a channel, both in 
bipartite and multiple access case one can consider a simplified BUC in general reasonings 



C(g)=pg+(l-p)UgU^. 



(1) 
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For two-access channels it holds g — g\ ® g 2 , where gi is the input of the i-th sender. 
2.2 Quantum error correction 

QECC is a subspace C of a larger Hilbert space %. Equivalently, a code is defined to be the 
projection Rq onto CCJf. One says that C is correctable if all states from this subspace 
g = RcgRc can be recovered after an action of a channel using some decoding operation T>, 
that is T> o C{g) = g. Such recovery operation exists if and only if 

R c A\A 3 R c = a i:j Rc (2) 
for some hermitian matrix — ctij. These conditions are due to Knill and Laflamme (KL) 

In the case of a larger number of senders we talk about local codes Ci, that is QECC for 
every sender. It is an immediate observation that KL conditions need only a little adjustment 
to serve for the case of MACs. Namely, we have (with the obvious notation): 
Observation 1 Local codes Ci are correctable for a MAC with Kraus operators {Ai} with k 
inputs if and only if 

Rci ® Rc 2 ® • • • ® R Ck A\AjR Cl ®Rc 2 ®---®Rc k = a ij R Cl ®R c ^®---®Rc k (3) 
for some hermitian matrix = otij. 

This is true since the set of product codes is a subset of the set of all codes. In further parts, 
we sometimes use the denotation R® R' (Rm <8> R'n) or S 53 S 1 for a code for a two-access 
channel and talk about M (g> N codes, where M, N denote dimensions of local codes. 

In case of many usages of a channel, A4 are tensor products of Kraus operators in KL 
conditions. In this paper, however, we concentrate on a single usage of a channel. For one 
use of a BUC, Eq. ([1]), KL conditions ((2]) reduce to the single condition (we write Rc shortly 
as R). 

RUR = XR, (4) 

which for MACs takes the form 



R®R'UR(g>R' = \R®R'. (5) 

It is useful to introduce the notion of the entropy of a QECC [29] . This entropy quantifies 
the number of ancillary qubits which are needed for the recovery procedure. For a BUC the 
entropy of a code is the von Neumann entropy of the matrix 



r - I P ±Xy/p(l-p) . . . 

I ±\y/p{T^pj 1-p I' 



By inspection one finds that S(C) = #(1/2(1 + yjl - 4(1 - p){\ - |A| 2 ))). It is equal to zero 
iff A = —1 or A = +1 and these values correspond to DFS with trivial (identity) recovery 
operation J 1 



a It is not clear to us whether it would make any sense to define entropy of a local code. 
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2.3 Higher rank numerical range approach to bipartite QEC 

Here we recall the notion of the higher rank numerical range f'and some of its implications 
in the area of QEC [161 117j . which are important for our current purposes. The approach 
we briefly describe below is the one we wish to further modify to be applicable in case of 
two-access channels. 

It is the form of KL conditions, which prompted the authors of Ref. [El [TT1 HU [20] to 
introduce the notion of higher rank numerical range (or the rank-fc numerical range) of an 
operator. For an operator A, it is defined to be the following set 

A fc {A) = { A g C : P k AP k = \P k } (7) 

with P k G Vk, where V k is the set of rank k projections. Elements of the set are sometimes 
called the compression values. 

Full characterization of the set A k (A) for hermitian A has been obtained. Namely, as- 
suming a,\ < a.2 < . . . < ajv is a spectrum of N x N hermitian A, it holds that A k (A) = 
(a k ,aN- k +i), which is (a) a true interval whenever a k < ajy-fc+i, (b) a singleton set if 
a k = ajv-fc+i, (c) an empty set in the remaining case. 

The set is also quite well understood for unitary operators [19]. We recall some results 
below. Let U be nondegenerate n x n unitary matrix with a spectrum spec(U) = 
corresponding to eigenvectors {|i>i)}™ =1 . Let A k (U) be the set of A such that for some k 
disjoint subsets d±, 6%, ■ ■ ■ , S k of {1,2, ... ,n} it holds A G conv{{zi,i G Sj}) for all i (conv 
stands for a convex hull). It was showrfthat A k (U) C A k (U). The proof of this fact is 
constructive in a sense that it gives explicitly the projection. Since 

with onj > and J2jeSi a n = 1 we can choose the code Pc to be Pq — J2i=i where 

= E v^j\ v j) ( 9 ) 

to obtain PcUPc — APc- One can take the subsets Si to represent triangles or, in a more 
restricted variant, sections. Both cases will be considered by us in Section [5] 

Relevance of the notion to the issue of construction of QECC is obvious [THl HZ] ■ 

3 Product higher rank numerical range and its basic properties 

Motivated by the form of KL conditions, Eq. ([3]) for multiple access channels we introduce 
the notion of the product higher rank numerical range. It is defined as follows 

Definition 1 The ki <8> <8 • • ■ product higher rank numerical range of an operator A is 
defined to be 

A kl ®k 2 ®-(A) = {A G C : R® R' ® • • • AR® Bf <g) • • • = Ai?® R' ® • • •} (10) 

6 It is the generalization of the notion of the numerical range of an operator X, which is defined to be the 
following set: A(X) = {A e C : (i/>|X|i/>} = A, <V#} = 1} 

c In fact for n > 3fc much stronger result was proved, namely A fc = A fc but to avoid technicalities we use here 
a weaker version, which is sufficient for our purposes. 
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for some R e Vk x , R' G T'fca , • ■ • 

In the above, we assume that it holds ki ^ 1 at least for a single index i . When all fcj's 
are equal to unity then one deals with the local (or product) numerical range [3 lj . that is a 

set 

Ai oc (A) = {A G C : {i/;\{cj>\A\^)\4>) for |$ G |0) G U b }. (11) 

If we choose all projections to be the same then we deal with the symmetric product higher 
rank numerical range A^™™^ (A). On the other hand, if only some of projections are the 
same we call the set locally symmetric product higher rank numerical range A k °^ s k v 2 ^ m . (A) . In 
the latter case, in a concrete case of interest, one needs to specify which projections are to be 
chosen the same. Obviously, the notion is meaningful only for operators acting on C ® C d . 
One has: 

Fact 2 It holds 

A m ® n (A) = |J A™ (1 ® UA1 ® C/t) = |J A%™-(V® lAV^ ® 1). (12) 

U V 

Proof: The defining equation (fTU)) can be rewritten as [R <g> (C/ti?C/)]A[i2 <g) (ZTtiJf/)] = 
\{R®WRU) for some unitary U. It is equivalent to {R®R){t®U)A(t®W){R®R) = \R®R 
from which the result follows. The second case is shown in a similar manner. ■ 
Naturally, one can rewrite the fact for the multipartite case. 

We propose to call the k\ ® k% (8 • • • product higher rank numerical range multipartite as 
opposed to the bipartite one, which corresponds to the k± <g> k,2 case. 

We also propose to use the dual set, the fci ® &2 ® • • ■ product codes set for A, defined as 
follows 

V kl ® k20 ...(A) = {ReV kl , R' E V k2 ,--- : R® R' ® ■ ■ ■ AR® R' ® ■ ■ ■ = \R® R' ® • • • }. (13) 

In a standard, i.e. non-product case, there is more than one projection (it may even be 
an infinite number, when the spectrum of an operator is degenerated) corresponding to the 
same compression value. This degeneracy may be removed in a product case, but we do not 
know to what extent this happens in a generic case. 

Although we will be interested mainly in BUC for which KL give the single equation ([5]) 
to solve, it is natural to introduce the notion of the joint product higher rank numerical range 
(just as it is defined in the stadard case [30]). Namely, one defines 

Definition 2 The k\ ® k% ® ■ ■ ■ joint product higher rank numerical range of operators Aj, 
i = 1,2, ... ,1 is defined to be the following set 

Afe!ig)& 3 ig)...(j4i, A%, . . . , Ai) = 

{(Ai,A 2 ,...,A/) G C 1 : Rig) R' ® ■ ■ ■ A,R® Rf ® ■ ■ ■ = XiR ® R' ® ■ ■ ■ } (14) 



for some R £ V kl , R' G V k2 , . 
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For further convenience, we also introduce the common higher rank numerical range which 
is defined as follows 

Definition 3 The ki ® &2 <8 • • ■ common product higher rank numerical range of operators 
Ai, i = 1,2, ... ,1 is defined to be the following set 

A fe°0fe 2 0...(^i^2,---,A/) = {A G C : R® R'® •■• A t R®R'® ••• = \R® R' ® ■•■} (15) 

for some R G Pki i R' G 7~ > k 2 >■•■ 

In the present paper, we will be interested solely in the cases when the product structure of 
the projectors corresponds to the tensor product structure of Hilbert space (which is usually 
uniquely determined by the problem under consideration). 

What differs the standard and the product higher rank numerical ranges most strikingly 
is the fact that while the former is determined solely by the eigenvalues of an operator, the 
latter would also be affected by the form of eigenvectors (in applications to construction of 
QECC it is enough to consider normal operators). This makes the product range difficult to 
determine even with the knowledge of the standard one. 

Product higher rank numerical range bears the following natural features |? 
Property 3 A m(g)n C A 

rnn 

Property 4 A m( g, n (A) can be empty even when A mn is non empty. 
Property 5 A m(gm C A loc 

Property 6 A m ® n (A) and V m ® n (A) are both compact sets 
Property 7 A (mi+m2)(g>n (A) C A mi(gin (A), i = 1,2. 
Property 8 A m ® n (A ® B) = A m (A) M A n (B). 

Property 9 A m ® n (A) = A ml $ n (U ® VAW ®V') for arbitrary unitary U and V . 

Property 10 A~(A,B) C A^ n (A, B) 

Property 11 A^™^ p (A) C A^^™'(A) C A m(g>ntg>p (A) 

Before we proceed, we need to recall the notion of the C-numerical range Wc(A) of an 
operator A, which goes as follows [32] : 

W C (A) = {trCtrfAU, U G It}. (16) 

When Us are taken to be product one deals with the local C-numerical range W^ C '(A) 33 . 
In cases considered in the present paper the set is a closed interval. With this notion in hand 
we can give the following bound on the product higher rank numerical range: 
Observation 12 It holds that 

A kl ® k20 ...(A) C w%g H , 9 ... (j^rA ■ ( 17 ) 

Proof: It is enough to take trace of both sides of the defining equation (fTU|) of the product 
range to conclude that every A belonging to A^®^®... (A) must also belong to W^ RI(g> ... {Ajk\ki ■ 



d In what follows wc use S to denote the Minkowski product of two sets on the complex plane, which is defined 
as follows: Z\ KI Zi = {z : z = z\Z2\ z\ £ Z\\ Z2 £ Z2}. 
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Naturally, we assume that W is local according to the same cut as the higher rank 
range is product. 

The observation will turn out to be very useful in bounding the product higher rank 
numerical range. Its application, however, will require a numerical optimization. 

We single out also some chosen properties of the common range, which will serve as a 
basis for one of the examples. 

Property 13 A% omm -(A, 1 - A) = {|} or 0. 

Property 14 Let Q — <S>Qi, with orthonormal basis {\i)} and rank qi projections 

Qi, act on d\®d 2 . Then following holds: (i) if I — d\ then Ad^kiQ) — A™ mm (Qi, Q 2 , ■ ■ ■ , Qdi)> 
(ii) if I < d\ then either Ad^kiQ) = or Ad^kiQ) = {0} with the latter holding if and only 
if0eA~{Q 1 ,Q 2 ,...,Q l ). 

Property 15 A£ omm -(Ai, A 2) ■ ■ ■ ,A K ) C A k &i A i) with E;^ = 1 - In particular, for 
K = 2, A—\A ll A 2 ) C A k (aA 1 + (1 - a)A 2 ). 

Proof of Property [13J Adding (i) RAR = \R and (ii) R(t-A)R = \R we obtain R = 2XR, 
thus A = 1/2 if equation (i) has a solution with this value, if it does not — Ak(A) is empty 
and so is A c ° mm -(A, 1 - A). ■ 

Proof of Property I14t (i) We look for A obeying 

R®R' (^>2\i){i\®Qi\ R®R' = \R®R' (18) 

for projections R and R' rank, respecitevly, d\ and k. Since Ejii NX* I = Idu ^ must be 
R = Id!, thus Eq. (Tr8"|) turns into 

Id! ® = Al rfl (g) i?', (19) 

from which immediately follows 

R'Q.R' = XR' Vj (20) 

concluding the proof of this part, (ii) replacing upper limit in the sum with I < d\ and 
considering diagonal terms we obtain 

R'QiR' = XR', i = l,2,...,l (21) 

and d\ — I equations = XR' . Thus A must be equal to zero and this value must be in the 
common range of Qi, i = 1,2, . . . ,1. The claim follows. ■ 

Proof of Property I15t Let K = 2 for simplicity. We have (i) RA\R — XR and (ii) 
RA 2 R = XR. We multiply the first equation by a, the second by 1 — a and add such 
equations to obtain (iii) R(aA\ + (1 — a)A 2 ) = XR. Thus each A which fulfills (i-ii) for some 
R, fulfills also (iii). The result then follows. ■ 

Our special interest in finding product higher rank numerical ranges of projections stems 
from the special type of noise we mainly focus on in the paper (see Eq. (|22[1 ). 
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4 Two access quantum communication: QECC for a BUC 

We now specify the type of noise we will further consider. We assume that U is hermitian. 
Hcrmiticity of a unitary matrix implies that it must be of the form 

U = P-Q (22) 

with some projections P, Q. 

4-1 Zero entropy codes (X = ±1) 

In Ref. [21] the following theorem was provec] 6 

Theorem 16 A M®N DFS for a channel Eq.^Tj) with the noise model Eq. \22\) exists if 
and only if at least one of the subspaces V = PH. and Q = QH is (d—M, d— N) -decomposable. 

It was also shown that it can be generalized to the generic case of U allowing for DFS, i.e., 
U = P + E k Qlh Qk with r(P ) > MN. 

In Ref. [21] we have applied the theorem to several concrete examples. We recall two of 
them (in both cases Q = span{|^>j)}) 

• for the following choice of 4>% no 2 ® 2 code exists 

|&> = ^(|11> + |22», 
|fc> = i(|10> + |21», 

• for the following choice the 2 ® 2 code does exist 

|^} = -L(|02} + |10», 
|0 2 ) = -L(|O1) + |2O)). 

4-2 Higher entropy codes(X ^ ±1) 

We now move to the case of higher entropy codes, that is for A = ±1 in Eq. ([5]). 
Let us start with some general remarks concerning the non product case, that is Eq. 
Assuming Eq. ([22]) and denoting 

7 =(l-A)/2 (23) 

we obtain 

RQR = jR, (24) 

or, equivalently, 

RPR = (1 - j)R (25) 



e See Ref. 1211 for explanation of the notion of the decomposability and the justifications of conclusions 
concerning examples. 
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It is almost evident that if Eq. (|24|) is to be fulfilled the ran l^of Q must be at least as large 
as the rank of R but it cannot be too large since Eq. (f25j) must also hold. The theory of the 
higher rank numerical range |16[ 117) makes this intuition strict and we have that < 7 < 1 
only when r(R) < r(Q) < n — r{R), where n is the dimension of the whole space. We will 
later prove a result of this kind for ranks of the reduced matrices of P and Q. It is clear that 
in such situations if 7 £ A(Q) then it is also the case for 1 — 7 for some other projection. 
In this way we have a dual pair of projections, say R and R. Notice that both codes have 
exactly the same entropy (see Eq. (j6|)). 

The main result of this section is the necessary condition for 7 to belong to the product 
higher rank numerical range of a projection. 

Theorem 17 Let Qi be a rank I projection. Assume Rm <8> R' n QiRm ® R'n — iRm ® R'n 
holds. Let further X\ > X2 > • • • be eigenvalues of trsQi, and yx > y% > . . . of tr aQi- Then 



M r(tr B Ql) 

^\xi-N~t\+ x i ^ MN V0- - 7)(1 + 37) + 1 ~ MN, 

i=l i=M+l 
N r(tr A Qi) 

^| yi -M 7 |+ x i ^ MNy/(\ - 7 )(1 + 37) + / - MN. 

i=l i=N+l 



(26) 



The proof of the theorem is based on several lemmae (some concerning also the non-product 
case), which we shall prove below. 

Lemma 1 Fix a number n. Let R = Yli=i l e i)( e i| ^ e a projection operator (the states 
\ei) are orthonormal) onto a subspace of an n dimensional space %. Denote with Br = 
{|ei), \e2), ■ • ■ , |e/c)} and B 1 - its orthonormal complement, so that span{Sfl iffiB 1 } = T~L. If Eq. 
\21$ holds with < 7 < 1 then Q must have the following form in Br® B 1 - , 



Q = 








^7(1- 7) 








S 



(27) 



where non zero blocks with off diagonal terms equal to zero are of size kxk and S is a projector 
of size (n — 2k) x (n — 2k) . 



Proof: If Eq. (|24[) holds then Q, which is hermitian, must necessarily have the form 



Q 



7lfc A 
At B 



(28) 



^In this context the rank of a matrix is understood in a standard sense. 
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with tk in the basis of the eigenvalues of R. Here A is a rectangular k x fc' matrix, while B 
denotes a positive semidefinite square matrix of size fc'. Since it is a projection it must be 
true that Q 2 = Q. For convenience we explicitly write the left hand side of this equation 



Q 



j 2 l k + AA* jA + AB 
7 At + BA* A^A + B 2 



Comparison of Eqs. (|28|) and (|29| gives us the set of conditions 

7 2 l fe +A4t = 7 l fc , 

A f A + B 2 = B, 
jA + AB = A. 

From the first one we obtain 

AA^=x 2 l k , x 2 =7-7 2 >0. 



(29) 

(30) 

(31) 
(32) 

(33) 



Denoting A = A/x we can rewrite this as A A* = Ifc. Thus AAJ must be rank k as r(lfe) = fc. 
Since for an arbitrary A it holds that r(AA^) — r(A) and r(A) < min{fc, fc'} we infer that 
fc' > fc as otherwise we would get a contradiction. Moreover A is an isometry so we have 



A = xA 



x^2\ei){vi\ 

i=l 



(34) 



with some orthonormal \vi) £ C fe . There are fc' such states in total but here only fc of them 
are important. We thus have 



AlA = x t Y l \v i ){v i \. 



(35) 



From Eq. (|3ip we conclude that [.A' -A, B] = since an operator commutes with its own 
function. The above equation implies that B has the following form 



B= ^A*|«i)(«j|J © B' 



(36) 



where B' lives on the subspace spanned by some orthonormal set The set of fc ar- 

bitrary non negative real numbers Xi needs to be determined. In the basis {|ei), \e%), . . . , |efe), \v\), . 
\vk)i \vk+\), . . . , \vk<)} the matrix Q reads 



7 


X 





7 


X 




X 


Ai 





X 


Afc 










B' 



(37) 
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Thus Q = Q' © B', where Q' has the following structure 











2 




■ : 





Since Q is a projection so are Q' and B' , which means that 



7 x 

x A fe 

7 x 

x Xi 



(38) 



must be a pure state. 



It can only happen if Aj = 1 — 7 . This concludes the proof as we take B 1 - = { | V\) , \v2) , • • - > | v ' k ) } 
and recall that x — y7(l — 7)- B 

Lemma 2 If RkQiRk = jRk, k < I < n ~ k holds for some < 7 < 1 then \\Qi — 7-Rfc||t r = 
V(!-7)(l + 37) + 



Proof: Due to Lemma [T] we have 





Qi - iPk 




(39) 



The matrix has a structure of a direct sum and we immediately obtain 
\\Ql - jRkWtr = k 



y/j(l- 7) 



+ l-k = ky/(l - 7 )(1 + 3 7 ) + I - k, 



vt(i - 7) 1-7 

where the last term is just the trace of S, which is rank I — k. I 

Lemma 3 Let the following hold for some 7 > and the rank I projection Qi 

Rm ® R'nQiRm ®R' n = jRm ® (40) 



TTiera: 



r(tr B Q z ) > M, r(tr A Qj) > N, 



(41) 



||tr B Q! - lNR M \\ tr < MNy/(l-y)(l + &y) + l " MJV, 
||tr A Qj ~ 7^'ivlltr < AfJVV(l- 7 )(l + 37) + Z - MTV. 



(42) 



Proof: We prove the result for one partial tracing as the second is analogous. For the proof 
of the first part we write for some arbitrary \<p) from TZm (a subspace with projection Rm) 



(ip\tr B Qi\ip) = J2 M(^I<WIV>;>+ (<P\{i>i\Qi\'P}\^i)- 



(43) 
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Recall now that Eq. (|4Q|) is equivalent to 

(<Pi\(lpj\Ql\<Ps)\lpm) = 8 is 8j m X (44) 

with orthonormal basis \<fii) and \tpi) for TZm and 1Z' N respectively. In virtue of this fact hrst 
sum is exactly NX. Second sum is non negative since Qi is positive semi-definite. It thus 
certainly holds that 

(ip\tv B Q l \i P }>NX>0, V MeRM . (45) 

Since matrix multiplication cannot increase rank, one has r(trsQ/) > t{Rm^bQiRm)- 
Thus it remains to show that Eq. (|4"5j) implies t{Rm^bQiRm) = t{Rm)- Take Rm = 
YhLi l^iXVil) tncn naturally Rm^bQiRm = Z)i=i 7« > °- It is now sufficient to 

put this into Eq. (|45|) to conclude that 7$ > for all z, which ends this part of the proof. 

The second assertion follows directly from Lemma [2] since for any G, H it holds that 
\\tr B G-tr B H\\ tr < \\G - H\\ tr [34]. ■ 

Proof of Theorem I17t The result follows directly from Lemma [3] and the fact that for 
Hcrmitian matrices A, B with eigenvalues a\ > ai > . . . > a n , bi > b% > . . . > b n respectively 
it holds that ||A-£|| tl > Yh=i l a * ~ bi \ E9- B 



Theorem [IT] can be also applied to — l«p ~ Qi giving the dual estimations 

Theorem 18 Let P&-i be a rank d — I projection. Assume Rm ® R'^PcP-iRm ® R'n 

(1—^/)Rm<EiR / n holds. Let further xi > X2 > • • • be eigenvalues oftiBPd 2 -h andyi > y% > • 

ofti A P d 2_i. Then 

M r(tr B P d 2_ l ) 



^|z;-iV(l-7)|+ ^ < MN V7(4 - 37) + d 2 - I - MN, (46) 

i=l i=M+l 
AT r(tr j4 P d2 _ ! ) 

\Vi - M(l - 7)1 + 51 il - MN y/l(* - 3 7) + d 2 -l- MN. 

i=l i=N+l 

By a direct inspection of the concrete case, we will show the following in the next section 



Theorem 19 Existence of a product code for the noise model U = P — Q for X = 1 — 27 
does not necessarily imply the existence of a product code for X = 2"f — 1 . In other words, 
there are cases when 7 € but 1 — 7 ^ ^2^2- 

This result provides sharp distinction between the standard and the product numerical range. 
Notice that Lemma Q] implies that eigenvectors of Q must have a special form. 

Corollary 1 If RQR = jR holds then there exists the basis in which k eigenvectors of 
the projector Q = J2i=i IV'iXV'il can be expressed as 



\<Pi) = V7l e i> + V /l -7l«i) 5 



(47) 
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where (e,|ej) = Sij, (vi\vj) = 5^, and (ei\vj) = 0. The code is then R = | ) (e^ | and its 
existence is equivalent to the existence of the code R — ^ 4 satisfying RQR = (l—~/)R. 

This implies that for any projection Q both 7 and 1 — 7 belong to Ak(Q). 

Thus if one wants to find a code for 7 7^ 0, 1 one needs to find a basis for which Eq. (|47|) 
holds. 

In Ref. [21) it was shown that in the case of qutrit inputs (d = 3) a zero entropy code 
is unique in a sense that there are no codes for the noise model Eq. (|22j) simultaneously for 
A = +1 and A = — 1. Higher dimensional codes were not considered there. With the above 
results in hands we can prove that this uniqness is stronger. Namely, we have 

Observation 20 If for d — 3 there exists a product 2® 2 decoherence free subspace then there 
exists no higher entropy code for this system. 

Proof: We will prove the result for a four dimensional projection operator Q4 since the result 
is proved in a similar manner for a five dimensional projection. We only show the part which 
has not proved in Ref. [3T]. Suppose equality S ® S'Q^S ® S' = S ® S' holds. This means 
that Q4 = S ® S' . The following equality R ® R'Q^R ® R' = 7.R ® R' with < 7 < 1 cannot 
then hold since this requires RSR = 7 and R'S'R' = 7' with 77' = 7. It is impossible since 
this implies that either 7 or 7' (or both) is not equal to one, which stays in contradiction to 
Lemma [T] (applying the theorem to S or S' we see that they must be at least four dimensional 
which contradicts assumption d = 3). 

Suppose now that R ® R'Q^R ® R' = holds. This implies that P5 := lg — Q4 = 
i? <8> i?' + with £ _L i? (g) R' . If there existed higher entropy code, it would hold that 

T <E) T'(R ® R' + |0 ®T' = 1 T®T' for some two dimensional projections T, T'. Hence 
T 1 ® T'|^)(^|T <E) T' must be a non zero vector if there is to be non zero 7. Otherwise the 
operators TRT and T'R'T' would have to be proportional to T and T" respectively which is 
impossible. This is a two qubit problem because of two dimensional projections so we rewrite 
it for clarity as E&F+ |*)(*| = 7I® 1. Thus G := 7I® 1 — i?®.F must be one dimensional. 
Let (ei,e2) and (/i,/2) be the eigenvalues of E 1 and i 7, respectively. The eigenvalues of G 
then read (7 — ei/i, 7 — ei/2, 7 — 62/1,7 — 62/2). Exactly three of these eigenvalues have to 
be equal to zero, which is impossible. I 

Finally, notice that the reasoning can be applied to an arbitrary noise model of the form 
U = P + e il3 Q. 

5 Applications 

In what follows, we will be mainly interested in finding bipartite product higher rank numerical 
ranges of operators acting on d ® d, one example will be devoted to the case of different local 
dimensions. 

5.1 Random swap 

Among all unitary matrices of the form (|2"2"|) . SWAP is the quantum gate probably most 
often used in the theory of quantum information. It can be written as SW APd^d = Psym. — 
Pasyrn. • — ^d: where Psym.-, Pa sym. are the projections onto symmetric and antisymmetric 
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subspaces respectively. Its action on pure states is given by Vrf 1 0) | V) = \ip)\4>). 

This kind of noise can be approached directly without resorting to the results presented 
above. We need to solve 

R ® R'V d R ®R' = \R®R' (48) 
with R = YliLi I fiii^Pi I an d R' — 52j=i I'^iXV'.jl- I* ls equivalent to 

(<Pi\{ipj\Vd\fk)\^i) =Mik8jl, (49) 

which in this case results in 

(<Pi\ipl)(<Pk\ipj)* = XS ik Sji. (50) 
As one can realize from the above A = is the only possible value. The corresponding 
eigenvectors of R and R' must obey ifi _L ipj , thus 

A 2 ® 2 (^)=0 (51) 

and 

Ajif»jv(Ki) = {0}, Af + A<d, d>4. (52) 

In the latter case with M = N = 2 we could take for example R = |0)(0| + |1)(1| and 
i?' = |2)(2| + |3)(3|. 

The above also means that we have: 

A£SOT) = 0. (53) 

Notice that the local numerical range (see Section [3]) provides a quick rough upper bound 
(0, 1) on the M ® N range. This is an example of application of the Property [7] 

Another bound narrowing previously mentioned one can be readily obtained from Fact 
Q21 which in this case does not involve optimization. Taking the trace of both sides of Eq. 
dH]) and exploiting the fact tvV d A (g> B = tr AB we obtain trRR' = MNX. Utilizing now 
the property (following from the Holder's inequality [36]) tr^4i3 < tiA holding for A > 
and I > B > and the fact that trAB > for positive semi-definite A and B, we obtain 
< MNX < min{Af, N}, which finally leads to 

< A<min{^,l}, (54) 

giving a significant improvement over the previous bound. 

Presented reasoning shows that in some cases direct exploitation of properties of U may 
prove very useful. 

5.2 Other kinds of noise 

5.2.1 The 3 ® 3 systems 
We start with d = 3 examples. 
Example 1. 

We will establish an outer bound for A2(g>2 of a projector Q(a) having the following eigenvectors 



^ 1 ) = (^|0) + VT r a|2))®|0), (55) 

|^> = |0)®(V5|1)+V^^|2)), (56) 

|V&> = |1> ® (v^|Q> + VT=^|2)), (57) 

|V> 4 > = (v^|l} + vT^|2»(8> |1) (58) 
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with < a < 1. From Corollary Q] we see that a £ A-2®2(Q(<x)) as all above vectors are of the 
required form y/a\ei) + \/l — a\vi). Notice that it also follows that a <E A2^ m '(Q(a)). 
Simple calculation yields 

/ 1 + a y/aQ - a) \ 

X = tr A Q(a) = tr s Q(a) = l + a y/a(l - a) (59) 

V y/a(l - a) y/a(l - a) 2 - 2a J 

and 

Y = tr A (1 9 - Q(a)) = tv B (1 9 - Q(a)) = 31 3 - X. (60) 

The eigenvalues are: X\ = 2, X2 = 1 + a, £3 = 1 — a and X\ = 2 + a, x 2 — 2 — a, and x 3 = 1. 
Setting M = N = 2 we get from Theorems [T7l and [TBI two inequalities on 7: 




|2 - 2 7 | + |1 + a - 2 7 | + 1 - a < 4^(1 - 7 )(1 + 3 7 ), (61) 
|2 + a-2(l- 7 )| + |2-a-2(l-7)| + l < 4^(4 - 3 7 ) + 1, (62) 

from which we get 

7eA 202 (Q(a))C^(l- A /l-^ a 2]^|i + 1 / 4 -3(^) |), (63) 

For each value of a we thus get an interval outerbounding One realizes that this 

example already provides the evidence in favor of Theorem II 91 

Using Observation 1121 this bound can be improved extending the region where a belongs 
to the range but 1 — a does not. As we have noted in the previous section, this usually 
requires optimization (in opposite to the application of Theorems IT7l i fTSl as above). We need 
to calculate 

A t (a) := iinintr[(C/(g)V)(Po«)^o)(C/ t «)^ t )Q(a)] (64) 



and 



A;(a) := - maxtr[([/ ® V)(P ® P )(U^ ® V^)Q{a% (65) 

4 u$gv 



where Po = |0)(0| + |1)(1|, which give rise, according to Observation [T2l to the bound 

A 282 (Q(a))C[A T (a);A 4 (a)]. (66) 

With this aim we parameterized U and V according to Ref. [37] and then, for different values 
of a £ (0, 1) with the step 8a = 0.01, computed A^ and Aj,. In this case the same bound is 
valid for A^™ ■ The results of this part of the paper are shown in Fig. QTsee also Section 
[7J. 

Example 2. 
Consider now the following operator [38] 

i<j 



Title 



A 2 ® 2 (Q(a)) 



1.0 




0.8 
















0.6 








* i a=0.7-> 




0.4 


















0.2 









0.2 0.4 0.6 0.8 1.0 



Fig. 1. The light shaded region represents the bound given by Eq. J63D in the whole range of a, 
the dark grey region is the bound given by Eq. Il66|l (the dots are for A-f- j_(af)). We have also put 
the section bounding the product range for an exemplary value a = 0.7. Additionally, wc plotted 
the lines a and 1 — a. While the former value always belongs to the product range (as we have 
discussed in the main text) there is a range of a, where the latter lies outside it with certainty as 
it is cut off by the bound, this cut off is represented by the dashed lines. 
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where 



= -jp= (ya\ij) + \/a~ T |ji)) , a + a 1 



d. 



We set d = 3. Then, explicitly, we have 



Q 4 =3 



/100010001\ 
OaOlOOOOO 
OOaOOOl 00 
010±00000 
1 1 1 
OOOOOaO 10 
OOlOOOiOO 
0000010±0 

V1000100 01/ 



with a — |(3 — \/5)- The reduced states are 



A,B 



tr 



g 4 = 



|(4±V5) 
I 







and 



Y AiB ee tr B)A P 5 = tr fl)A (l9 ~ Qi) = 31 3 - 



H 4 TV5) 



i (5 t V5) 



i(5±V5) 

from which the inequalities on 7 (Theorem [T7] and its dual Theorem ITS)) follow 



i(4 + 75)-2 7 
^ (5 + Vg) - 2(1 - 7) 



2(1-7) 



+ I (4 - Vs") < V(l- 7 )(l + 3 7 ), 



1 



5- V5) < 4\/7(4-37) + l. 



This eventually leads to the following outer bound 

4-V13 3 + V33 



A 



2® 2 



c 



[0.0657; 0.972]. 



6 ' 9 

At the same time, application of Fact [T2l gives us better estimates: 

A 2 ® 2 C [0.1788; 0.7378]. 



(69) 



(70) 



(71) 



(72) 
(73) 

(74) 
(75) 



We have also used Observation [12] to estimate bounds on symmetric product range and 
we have obtained that 

A^ m ' C [0.4166; 0.5556]. (76) 

We were not able to verify whether any of the values inside these intervals actually belong 
to the (symmetric) product higher rank numerical range of Q. 
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5.2.2 The 4 ® 4 systems 

Now we give an example of the projection Q for which the numerical range A2&2 can be found. 
New analytical technique will be introduced with this aim. Let the eigenvectors of Q(a) be 

1^) = v^|00) + VT - a|22), (77) 
IV*) = v^|01> + vT^|23), (78) 



|V' 3 ) = v^|10) + Vl^|32), (79) 



|V»4> = v^lH) + VT~c^\33). (80) 

Theorem [1] leads to the conclusion that both a and 1 — a belong to A2®2(Q(aO). Assume 
1 — a > a. Consider now projectors i? and R' spanned respectively by 

|0 1 )= A /W|O) + V^|2), (81) 

|«fe) = Vl - + v^|3) (82) 

and 

l^i) = |2>, (83) 

|^ 2 ) = |3) (84) 

with j3 € (0,1). One can easily check that with such choice of projectors we can get 7 = 
/3(1 — a) and thus varying /3 over the whole range in turn that 7 G (0, 1 — 7). We will now 
show that this is best what one can do and this interval represents A2®2(Q(oO)- 

With this aim recall once again that if 7 € A202 then there exist projectors R = 
l&Wil + l<fe)<<fe| and i?' = |V>i)<Vi| + 1^2)^21 such that (^(^IQMHi) = jSikSjh 
Naturally (<fii\4>j) — Sij and = Sij. Let the vectors be decomposed as follows 

\<t> 1 ) = \e,e'), e,e'eC 2 , (85) 

|0 2 ) -I/,/'), /,/eC 2 , (86) 

\i, 1 ) = \E,E'}, E,E'eC 2 , (87) 

|V 2 ) = |F,F'), f,f'£C 2 , (88) 

where |i, j) is understood as the direct sum |i)©|j), i.e., we use the isomorphism C 4 = C 2 ©C 2 . 
Let now p,p' , q, q' , P, P' , Q, Q' be variables in C 2 , which can take values as follows p,q = e, /; 
p\ q> = e ', /'; P,Q = E, F; P' , Q' = £7', F' in combinations consistent with Eq. ([55]). After a 
direct calculation we arrive at 



(p,p'\(pp'\Q4hq')\Q,Q') = (v^pK-PI + VT^WQ (V5|«)|Q) + VT^WM)) , 

(89) 

which must be zero whenever (p,p') ^ {q,q') or (P,P') ^ (Q,Q') and 7 otherwise. Recalling 
what values can be taken by the variables and orthogonality conditions we thus conclude 
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Lemma 4 The number^/ belongs to A-2®2{Q(oi)) if and only if there exists a set of four vectors 
from C 2 <g> C 2 



\ X i)=V^\e)\E) + VT :r ^\e')\E>), 
\ X 2) = Me}\ F ) + VT^\e}\F'), 
\ X3 ) = Mf)\E) +VT^\f')\E'), 

\ Xi ) = Mf)\F) + VT^\f')\F>) 



such that 



\\\E)\\ 2 + 
\\\F)\\ 2 + 
(e\f) 
(E\F) - 



\\E' 

\\F',,, 
(e'\f) 
(E'\F' 



1, 

1, 
1, 
1, 
0, 
= 



and 



(Xi\Xi) = lk 



This observation is crucial since now we can bound 7. Writing 
have for example for \xi) 



Vi = Vbaba) = llxill = \\V$e)\E) + VT^\e')\E' 

^ ||e|| • + VT^ He'll • ||£'|| < VT^ (||e 



< 



E\ 



\E\\,\\E'\\] < VT^ 



VT~a vj|e| 
which finally leads to 



+ l|e'HVl|£|l 2 + l|£'ll 2 = VI^, 



A 202 (Q(a)) = (0,1 -a). 



as 



(90) 
(91) 
(92) 
(93) 
(94) 



(95) 
(96) 
(97) 
(98) 
(99) 
(100) 



(101) 

|i| I for clarity, we 



\E\ 



\E'\\ 
\\E'\ 



(102) 



(103) 



Let us just explain all the important steps in the estimation for 7: (i) the first inequality 
follows from the triangle inequality, (ii) the second inequality uses a < 1 — a, (iii) the third 
is the Cauchy-Schwarz inequality, (iv) the last equality makes use of normalization. 

Proceeding in a similar manner in the case of a > 1 — a we find that in this region the 
product numerical range is A 2( g, 2 (Q(a)) = (0,a). 

Interestingly, inclusion of sets from Observation [T^] becomes an equality in this case so 
the same region can be found by numerical optimization of the bounds on the local R £§> R'— 
numerical range of Q(a) (see previous examples). 

It turns out that in this case 



A5g3 m -(Q(a))=A 282 (Q(o)) > 



(104) 
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what we prove below. The set equal to (a, 1 — a) if < a < h and (1 — a, a) otherwise can 
be shown to belong to A^g^" 1 ' (Q(a)) simply by taking ipi = 4>i in Eq. ()83[) and leaving <f>j 
unchanged. The remaining region can be shown to be achievable by setting 

= = v / W|0)+iV^|2), 
|0 2 ) = |^ 2 ) = v / W|l)+iV^|3). 

The results of this section are summarized in Fig. [2j in which we have also plotted the 
bound stemming from Theorems 1171 and 1181 

A 2 ®2(Q(<*)) 




0.2 0.4 0.6 0.8 1.0 



Fig. 2. The shaded region is the 2 ® 2 range of Q(ol) given by Eqs I|77I80|I . The hatched region is 
the bound from Theorems 1171 and 1181 



5.2.3 The 2 ® 4 systems 

Consider the following projector acting on C 2 (g) C 4 

Q( 7 ) = |0)(0| ® Qi + |1)(1| <g> Q 2 ( 7 ) (105) 

with 

Qi = |0)(0| + |l)<l|, (106) 
Q 2 (7) = |r7i)<?7i| + I%}<%|, (107) 



m = v7io) + yr^i2), 

l?72)=V7|l> + v / l — 7|3». 



(108) 
(109) 
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From Property[T4]it follows that we must find A2° mm '(Qi, Q2(l))- Taking a = 1/2 in Property 
13 we obtain A| omm -(Qi, Q 2 {l)) Q A 2 (Z( 7 )) with 



Z( 7 ) = 



1 + 7 








V 



^7(1-7) 

1+7 ^7(1-7) 

>/7(l - 7) 1-7 



\ 











Vt(i - 7) 







(110) 



1-7 J 



By inspection, we find the eigenvalues of ^(7) to be | (1— y/y), ^(1 — ^/7), 5(1+^/7), 5(1 + ^/7)- 
From the theory of higher rank numerical range (see Section [2. 3[) it immediately follows that 




Fig. 3. The shaded region is the product higher rank numerical range of Q("f) given by Eq. Q105I . 
The hatched region in which it is included is the bound stemming from the application of Theorems 
IT7landlT8l 



^(1-V7)^(1 + V7) 



(111) 



It remains to show that the above represents the sought common numerical range (and in 
consequence the product numerical range), that is all A G A 2 (Z(7)) are achievable in the 
sense of fulfilling the following equations: (i) RQ\R = XR and (ii) RQ2('j)R = XR for some 
R. Indeed, take R to be projecting on the subspace spanned by the following two vectors: 



|6) =V\\0)+c^VT^ A|2), 
16) = VA|l)+c^yi^A|3). 



(112) 



Then, by varying j3 in the range [0,7r], we can get any A from the desired interval. 
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We have plotted the range in Fig. |3] along with the bound stemming from Theorems 1171 
and[T8l In Fig. 0]we show application of Property [TBI for values a — 1/2, 1/3, 1/4, 1/5, 1/10. 

We have also verified by optimization that Observation [12] gives exactly (|111[) as the outer 
bound. 

A 2 ®2(Q(r)) 




0.2 0.4 0.6 0.8 1.0 



Fig. 4. The inner shaded region corresponds to a = 1/2 in Property 1151 and represents the true 
2 ® 2 rank of Q(-f) given by Eq. (11051) . The region with dashed lines corresponds to a = 1/3, solid 
lines — a = 1/5, dot dashed lines — a = 1/10. In the limit a — > we will obviously get [0, 1]. 

6 Reverse problem — a toy model 

We conclude with an example of a different type of noise. Unitary U which we will consider 
have a regularly distributed spectrum and can hardly be called generic, nevertheless, as it 
can be treated analytically, we believe that the examples may be an important toy model for 
future work on product numerical range. We concentrate on highest entropy codes, which 
correspond to A = 0, but the observation can be generalized to other compression values. Our 
examples may be considered as an illustration to the problem inverse to the considered so far. 
That is 

Problem. Given a product code R® R' construct a class of unitaries U for which 

R®R'UR®R' = 0. (113) 

We will focus on 2 (g) 2 codes and consider d = 4 case. Using the terminology introduced in 
Section [3] this problem can be rephrased as follows: given an element of a product codes set 
of some unitary operator U corresponding to the compression value A = 0, find an exemplary 
form of U. Not surprisingly homogenous equation on U is very different from homogenous 
equations on projectors considered in |21j so we need a totally different approach. 
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In our case we have P = R® R' , i.e., the sum iV'iKV'il lS a product projector (J 

is a set of indices). It is worth stressing that do not have to be product itself — they 
only need to span four dimensional product subspace. We now need to construct spectrum 
of U = Zi\vi)(vi \ so that A = is a compression value. This can be easily done if we recall 
the well known identity 

1+oj + oj 2 = 0; U3 = e^\ (114) 

We choose 

spec([/) = = 

{e^ 1 , oje^ 1 , w 2 e 1?1 , f#* , we' 6 , w 2 e 1?2 , e 1 ^ , we' 6 , w 2 e ifo , e'« 4 , ue iU , wV« 4 , e i£s , e lfo , e* 7 , e i?s } 

with arbitrary £j 7^ (at least for i,j = 1,2,3,4). As it was mentioned earlier, we can 
construct triangles from the eigenvalues, that is we can now take S m — {e^ m , we'^ m , w 2 e l£m } 
and due to Eq. (|114l) this means that indeed A = is a compression value since we can set 
the numbers a's to be all equal to 1/3 (see Eq. ((HJ). It remains now to properly choose 
eigenvectors. This construction uses the same arithmetic identity. We take \vi) to be such 
that all terms in superposition besides the first one cancel due to Eq. (|114j) when added with 
coefficients l/y/3 stemming from Eq. ([9]). The states remaining after the summation should 
sum to a product projector. We have explicitly in a closed form 

i = 1,4,7,10 (115) 

1+2) / 

with orthonormal states \ui) for z = 1,4, 7, 10 such that 

\ui){ui\=R®R'. (116) 

i=l,4,7,10 

We could, for example, take two-qubit Bell states embedded in a two ququart space, i.e., 

K, 4 ) = ^(|oo)±|ii)), 

|n 7 , 10 ) = -^=(|01>±|10». (117) 

The rest of eigenvectors may be chosen arbitrarily with the only restriction that all \vi) are 
orthonormal. 

Let us summarize for clarity all the elements. We have 

2 1 

V-^ +fe = 0, fe= M,7, 10, (118) 
*=o 6 

which ensures that A = belongs to numerical range. We then assume that Eq. (|116p holds 
so we can define 

I 2 

:= 75 J2 l u ™+ fc ) = fc = 1, 4 ,7, 10, (119) 

v m— 
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which properly sums to a product code. 

We also propose an alternative construction of eigenvectors. The following set of vectors 
\vk) allows us to construct a code 

V2e ia * \ / \ Uk ) \ 

_^ e ia fc ^ e i Qfc | Ufc+i) fc= 1,4,7,10 (120) 

_^ e i Qfc _\^ e i« fc J \ \u k+2 ) J 

with some arbitrary phases ctj and \uk) defined as previously. 

One can check that in the considered case of highly regular spectrum requirement about 
non degeneracy of spectrum can be relaxed and the phases for k = 5, 6, 7, 8 can be com- 
pletely arbitrary. 

Consider now a similar example in d = 3 in which, instead of triangles, we will consider 
construction of sections crossing in the z = point so that it belongs to the numerical range. 

We have U — J2j=i e ' aj with ct; > ctj for i > j, and a\ — as — 7r, a2 = a§ — 7r, 

«3 = «7 — 7r with an additional constraint that € conw(e 1Q4 , e 1Q5 , e 1QQ ). Fig. [6] shows an 
exemplary distribution of eigenvalues. 

In general: 

/ 1 \ / 1 \ 

U B = (l-\f 3 )(M)®(l-\w 3 )(w 3 \) (121) 



/I 














R®R' = U A 
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j U\®U B 1 
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with some unitary Ua and Ub and corresponding vectors I/3) and \w$). 
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Fig. 6. Exemplary distribution of eigenvalues on unit circle. We have assumed, as further in the 
main text, that 04 = 09 — 7T. Moreover, ce± = c*5 — n, a.2 = ctf> — it, 03 = 07 — n. The arrows 
represent the summation, Eq. I|123|l . giving rise to the compression value 0. 



Set now the Fourier matrix F3 in place of the local unitaries Ua and Ub- 



U a = U b = F 3 




uj = e 3 



(122) 



Then I/O = (l,w,w 2 ) T , |/ 2 ) = (lo 2 ,1,lo) t , |/ 3 ) = (c^ 2 , 1) T , \ m ) =\ fi ). We now define 
properly eigenstates of U, which could lead to a product code P = X)i=i iV'iXV'il = R® R' , 
Hi) = l/i> <8 |^ 2 > = ® |/ 2 >, |V>3> = I/2) ® |/l>, 1^4) = I/2) ® I/2), with a simplifying 
assumption = ag — ir. By choosing our eigenvalues so that 



1 



; z j + o*j'+5 = 0- .7 = 1,2,3,4 



(123) 



we ensured that we can set 



\1>i 



V2 



(|^) + |0 i+ 6», j = 1,2,3,4. 



(124) 
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It remains to choose \<pi) so that \ipi) are product. For example, one choice could be 



|^l,6> = i(|/l>C 


5 |/i) ± |/ 3 ) $ 


M/2)), 


(125) 


l^ 2i7 ) = -L(l/ 1 )c 


3 I/2) ± I/2) $ 


5 l/a)), 


(126) 


|03, 8 ) = ^(|/a) c 


^ l/i) ± l/i) $ 


9 l/3», 


(127) 


|04,9) = ^(|/a> 


8 I/2) ± l/a) C 


3 l/x»- 


(128) 


Our construction is in fact such that belongs to the symmetric product higher 


rank 



numerical range. 

7 Discussion and conclusions 

Motivated by the form of Knill-Laflamme conditions for multiple access channels, we have 
introduced the notion of the product higher rank numerical range as a tool helpful in con- 
structing quantum error correction codes for such type of quantum channels. Several useful 
extensions of it, namely: the symmetric product range and the common product range have 
also been discussed. Techniques for bounding the product range and some analytical tech- 
niques for findings ones in some cases have been introduced. We have applied our findings to 
construction of error correction codes for a class of two-access biunitary quantum channel. 
The reverse problem of finding the noise model for a given product error correction code has 
also been considered. 

Concluding, we state some open problems related to the subject. Among them, determi- 
nation of the shape of the product higher rank numerical range for different types of operators 
seems to be one of the most important ones. In particular, for the case of d = 3, it should be 
verified whether the 2® 2 product range of a projection operator is at most a one-element set. 
If this is indeed the case, the uniqueness (Observation [20)) is a general feature in this setting. 
Also the issue of simple-connectivity of the set Ak ll g ) k 2 g ) ...(A) in arbitrary dimensions is of 
a particular interest. From the point of view of quantum error correction special attention 
should be devoted to normal operators. 

We must note that in principle there exists a method for judging whether a given value 
A belongs to the product range of a given operator (and determining the corresponding pro- 
jection if it does belong) . It goes as follows (for simplicity we focus on the bipartite range) . 
We set Pq — |0)(0| + |1)(1| and parametrize unitary operators U and V. Then we compute 
D := min^v \\Po® PqU ®V AW Pq®Pq — APo<g>Po|| for some appropriately chosen norm 
|| • || [5S]. Clearly, if D cannot be made smaller than the prescribed accuracy then A does 
not belong to the range. However, it is obvious that we need something better for practical 
purposes. There is thus a pressing need for designing an efficient procedure for finding explicit 
forms of product projections from the product codes set. 
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